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Abstract. The chiral Luttinger liquid theory for fractional quantum Hall edge transport 
predicts universal power-law behavior in the current-voltage (I-V) characteristics for electrons 
tunneling into the edge. However, it has not been unambiguously observed in experiments 
in two-dimensional electron gases based on GaAs/GaAlAs heterostructures or quantum wells. 
One plausible cause is the fractional quantum Hall edge reconstruction, which introduces non- 
chiral edge modes. The coupling between counterpropagating edge modes can modify the 
exponent of the I-V characteristics. By comparing the i/ = 1/3 fractional quantum Hall 
states in modulation-doped semiconductor devices and in graphene devices, we show that the 
graphene-based systems have an experimental accessible parameter region to avoid the edge 
reconstruction, which is suitable for the exploration of the universal edge tunneling exponent 
predicted by the chiral Luttinger liquid theory. 



1. Introduction 

The fractional quantum Hall (FQH) effect has been observed in GaAs-based semiconductor 
heterostructures and quantum wells [1] for almost 30 years. The fractional quantum Hall states 
have a lasting fascination due to their nontrivial topological properties, which lie beyond the 
paradigm of the Fermi-liquid theory, and the potential applications in topological quantum 
computation [2j [3j Hj . Probing excitations at the edge of a FQH droplet is the most accessible 
approach for the detection of the topological properties in the bulk thanks to the edge-bulk 
correspondence in the topological system. The edge excitations are gapless in the thermodynamic 
limit, contrast to gapped quasihole/quasiparticle excitations in the bulk. Owing to the existence 
of edge states, current exists between two contacts connected by an edge channel as electrons 
can be injected into or removed from the FQH edge without costing energy. The standard 
theory for the FQH edge physics is the chiral Luttinger liquid (CLL) theory [5j. The theory 
predicts that a FQH droplet exhibits power-law bahavior in the I-V characteristics (/ oc V a ) 
when electrons are tunneling through a barrier into the edge, e.g., from a three-dimensional 
Fermi liquid. For the celebrated Laughlin state at a filling factor v = 1/3 in the lowest 
Landau level, the CLL theory predicts a tunneling exponent a = l/^ = 3[5]. However, the 
exponent a measured in experiments are sample dependent with a numerical value smaller 
than 3, which are nonuniversal [6j EJ [9] . One of the possible causes of this discrepancy is 
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Figure 1. Simplified models to mimic 
the experimental setups in (a) GaAs- 
based and (b) graphene-based systems. 
The GaAs model consists of an electron 
layer and a homogeneous impurity 
charge layer, which confines electrons. 
The graphene model consists of layered 
Dirac electrons with an open edge and 
their image charge due to metallic gate; 
both confines electrons. 



existence of counterpropagating edge modes, which result from edge reconstruction [§1 |1U| [TT]. 
The edge reconstruction in the electron density profile manifests the competition between 
electron-electron interaction, which tends to spread out the electron density, and the electrostatic 
confinement at the edge, which holds the electrons in the interior of the sample. Numerical 
calculations \V2\ [P5] indicate that the edge reconstruction in FQH states is unavoidable in the 
state-of-art GaAs-based samples based on modulation doping, which greatly enhances sample 
mobility but unintentionally introduces a large region of fringe field at the edge that favors edge 
electron redistribution. This effect hinders the observation of the universal behavior in edge 
transport. 

In recent years new two-dimensional materials, such as graphene, emerged. The FQH effect 
has already been observed in graphene |14 |I15[[16] . However, integer and fractional quantum Hall 
effects in graphene differ from those in normal semiconductor-based devices, because electrons in 
graphene obey the relativistic Dirac equation, instead of the nonrelativistic Schrodinger equation 
often encountered in solid state materials. Here we compare the FQH states in both GaAs-based 
samples and in graphene samples with the emphasis on the difference of the competition at 
the edge. We find that there is an experimental accessible parameter region in graphene-based 
samples, in which edge reconstruction can be avoided and, therefore, universal FQH edge physics 
can be observed. The universality is reflected in the calculation of the equal-time edge Green's 
function, which is related to the I-V characteristics in edge tunneling. 

2. FQH states in GaAs-based systems 

An experimental sample based on a modulation-doped GaAs/GaAlAs heterostructure |17| can 
be simplified to a two-dimensional electron layer, which is located at the interface between 
GaAs and GaAlAs, and a homogeneous positive background charge layer, which originates from 
the remote dopants. This is illustrated in Fig. Q] (a). The confinement potential for electrons 
due to background charge competes with the Coulomb interaction between electrons. This 
competition is the driving force for edge reconstruction. The density of the background charge 
a is determined by the electron filling factor v due to charge neutrality and, therefore, the 
background confinement potential is a single-body potential in the FQH Hamiltonian. The 
complete Hamiltonian can be written as 



where the c+ is the electron creation operator for the lowest LL single electron state with an 
angular momentum m. 
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Figure 2. (a) The lowest energy 
in each total angular momentum 
subspace for 8 electrons in 24 
orbitals with different ds in the 
GaAs model. The global ground 
state is labeled by an arrow, 
(b) The corresponding angular 
momentum of the global ground 
state as a function of d. 
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Figure 3. The low-energy spectrum for 8 electrons in 26 orbitals with (a) d = 0.5 Ib, (b) 
d = 1.0 Ib, (c) d = 1.4 and (d) d = 2.0 Ib in the GaAs model. Edge reconstruction, signaled 
by the ground state angular momentum change, occurs around d c ~ 1.4 Ib- 



represents the electron-electron interaction and 

U m = ea [ [ d 2 nd 2 r2 ___L = ==|^ m ( ri )|2 (3) 

J Jr 2 <R V" + Fl _ r 2\ 

the confinement potential. The distance d between the electron layer and the background charge 
layer is the parameter that controls the relative strength of the confinement potential to the 
electron-electron interaction. When d is small, the confinement is strong and electrons tend to 
stay in the region defined by photo lithography process (r < R for the disk sample we consider 
here). However, the confinement becomes weak when d is large and thus the electrons near the 
edge tend to spread out, leading to edge reconstruction. 

Because of the rotational symmetry in the geometry and the gauge we choose, the total 
angular momentum Mtot is a good quantum number and thus we can diagonalize the Hamiltonian 
in each angular momentum subspace. In Fig. [2(a), we plot the lowest energy in each of the total 
angular momentum subspace with different ds for 8 electrons in 24 orbitals. As discussed in 
Ref. [10\ 112] . when d is smaller than the critical value d c = 1.5 =t 0.1 Ib, the global ground state 
has the same total angular momentum Mtot = Mq = 3N(N — 1) /2 as that of the Laughlin state. 
When d > d c , edge electrons spread out under the weak confinement and, therefore, the total 
angular momentum of the global ground state increases as shown in Fig. [2(b). To further reveal 
how the edge reconstruction happens with increasing d, we plot the low-lying energy spectrum 
in Fig. [3] for 8 electrons in 26 orbitals with different ds. The larger number (26 rather than 
24 or 22) of orbitals is chosen to avoid unnecessary mixing of the bulk and edge modes in a 
finite system. The chiral bosonic edge mode (dashed line), which describes the edge electron 
density deformation, is linear for small AM = Mtot ~ Mo but develop downward curvature as 
AM increases |12j . As d increases, the bosonic mode bends down further and eventually the 
total ground state momentum jumps from Mtot = 84 to Mtot = 91 above the critical value d c ; 




Figure 4. The electron density profile for the ground state at (a) d = 1.0 Ib, (b) d = 1.4 Ib and 

(c) (f = 5.0 Ib in the GaAs model. The dashed line is the density profile for the Laughlin state. 

(d) The edge Green's function before and after edge reconstruction; the slopes in the log-log 
plot yield a = 3.2 and 2.1, respectively. 



hence edge reconstruction occurs. Fig. H^a)-(c) shows the evolution of the electron density profile 
for the ground state as d increases. For large enough d, the electron density profile develops 
an almost detached edge piece, as shown in Fig. Sfc), suggesting the reconstructed edge can 
accommodate three counterpropagating edge modes. The edge transport thus depends on the 
coupling of the edge modes, hence no universality can be expected except in the strong coupling 
limit. 

To explore the universality of the FQH edge states, we calculate the equal-time edge Green's 
function, 

? w^mpm^ (4) 
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where \E r+ (r) [or Vl/(r)] creates (annihilates) one electron at fat the edge of the FQH droplet. 
In the disk geometry with |r| = \r'\ = R, the distance \r — r'\ = 2i?sin(#/2) where 9 is the angle 
between r and r 1 . The CLL theory predicts that in the large distance limit, the edge Green's 
function has an asymptotic behavior 

\G{r-?)\ ~ \r- r'\~ a oc | sin(0/2)| _a . (5) 

For the v = 1/3 Laughlin state, the CLL theory predicts an exponent a = 3. Fig. [U[d) plots 
the edge Green's function for 8 electrons in 24 orbitals with Coulomb interaction for two ds 
on the opposite sides of the edge reconstruction transition. For comparision, the edge Green's 
function for the variational Laughlin wavefunction is also presented. We find an exponent 
a = 3.2 ± 0.2 [18] at d = 1.0 Ib < d c , which is consistent with the CLL theory prediction. The 
exponent for the variational Laughlin wavefunction is the same as that of the Coulomb ground 
state before reconstruction. However, when d = 2.0 Ib > d c , the exponent a = 2.1 is noticeably 
smaller than 3. This suggests that the long-range Coulomb interaction does not directly lead 
to nonuniversal edge behavior, but the consequence of its competition with the confinement 
potential, i.e., edge reconstruction, does. In a typical GaAs-based sample, d is of order 1000 A, 
i.e., about 10 Ib [17] . While the numerical calculation for such a large d is not practical due to 
the system-size constraint, we speculate that multiple edge reconstructions can occur, leading to 
more than three counterpropagating edge modes. Therefore, such a device is always in the edge 
reconstructed regime, which does not guarantee the observation of a universal edge tunneling 
exponent. 



Figure 5. (a) The single- 
particle lowest Landau level en- 
ergy in the graphene model with 
an open edge clS Si function of 
the angular momentum in a disk 
sample with 60 orbitals. (b) The 
angular momentum of the many- 
body ground state clS Si function 
of d for 8 electrons in 32 orbitals. 

3. FQH states in graphene-base systems 

The experimental setup for graphene samples is illustrated in Fig. [I^b) (and Ref. [21] ). The 
graphene layer is situated at a distance d from a metallic gate (boron-nitride) |16j . Therefore, 
electrons have their image charge at a distance 2d from the graphene layer. This modifies the first 
term in the Hamiltonian in Eq. (pQ), which now contains both the electron-electron interaction 
and the electron-image charge interaction: 
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The combined interaction behaves like dipole-dipole interaction with 1/r 3 dependence at short 
distances and thus the interaction pseudopotentials (V m s) [H] are dominated by those with 
small ms. The single-body term in Eq. ([I]) comes from the edge confinement potential in the 
graphene due to the open boundary. The solution of the Dirac equation for a semi-infinite 
graphene with an open edge at x = satisfies equation D^— \^2x c ) = [20] where D^{x) is the 
parabolic cylinder function and x c = Mb, A* = e 2 /2; k is the wave vector parallel to the edge 
and e the energy in units oihvp/lB- In the disk geometry with N or i, orbitals, the edge is located 
at r c = yJ2N or blB, thus the solution satisfies [— y/2(x c — Vc/Ib)] = 0. Solutions are allowed 
for u = 0, 1, 2, ■ ■ ■ in the bulk, hence the energy for the nth LL is e n = y2n. As plotted in 
Fig. E(a), the eigenenergy of the lowest Landau level increases rapidly as approaching the open 
edge, which keeps electrons from spreading out. Since in the disk geometry the mth orbital is 
located at y2mlB along the radial direction, the Landau orbitals closer to the edge are denser 
than those closer to the center. The edge potential influences a large portion of the orbitals 
in the momentum space. Therefore, in a finite-size calculation we need to include significantly 
more orbitals than needed for the Laughlin state to prevent the edge potential from destroying 
the bulk properties. 

In Fig. EJb) we plot the ground-state angular momentum for 8 electrons in 32 orbitals as a 
function of d. Note that for small enough d the ground state can have an angular momentum 
smaller than Mq = 84, in sharp contrast to the GaAs model in Fig. EJ^b), in which the Laughlin 
phase remains stable down to d = 0. This is because that the interaction is too weak to compete 
with the edge potential when d is sufficiently small, thus electrons are repelled from the edge and 
the ground-state angular momentum decreases. Strictly speaking, a fair comparison requires us 
to use a sharp cutoff in real space in the GaAs model as well, such that the single-particle energy 
bends up and one needs to introduce Landau level mixing in a more complete treatment (see 
Ref. [12J for more detail), which is irrelevant to our main results so we skip here. Essential 
to our discussion, the Laughlin phase is nearly system size independent when N or b is large as 
analyzed in Ref. [21]. As pointed out there, we believe that one has a good chance of realizing 
the Laughlin phase in graphene without reconstructed edge for 0.5 < d/ls < 1.5 [21]. The 
upper bound of the window in which the unreconstructed Laughlin phase can be realized in the 
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Figure 6. (a)-(c) The low- lying energy spectrum for 8 electrons in 32 orbitals in the graphene 
model with various ds. The corresponding Laughlin state has a Mtot = 84, the same as in (b). 
(d) The exponent a in the edge Greens function of the ground states for various ds. 

graphene model agrees well with that in the GaAs model |10| , I12 |, [T3] . However, in the graphene 
model, e.g., on boron nitride substrate between the graphene layer and the gate, the parameter 
d can be chosen as small as a few nanometers or a small fraction of Ib |16| . thus the universal 
exponent in tunneling experiments may be easily accessible in graphene. 

In Fig. E|a)-(c) we plot the low-lying energy spectra for d = 0.7, 1.1, and 1.7 Ib, whose 
ground-state angular momenta are 76, 84 (Laughlin), and 91, respectively. In these cases the 
edge mode is not well separated from other states in the spectrum (similar to Fig. 8 in Ref. [12]). 
To demonstrate the universal tunneling exponent, we compute the edge Green's function for 8 
electrons in 32 orbitals, as shown in Fig. [6^d). We find that a for the three ground states are 
quite different. They are 1.6 and 2.3 for d = 0.7 Ib and d = 1.7 Ib, respectively; neither is in 
the Laughlin phase. In these casses, the Laughlin state still suffers from the edge reconstruction 
instability, which results in charge accumulation similar to that studied earlier based on Hartree- 
Fock types of calculations |22| ,l23[ [2"4" l l25]. When d = 0.9, 1.1 and IAIb within the Laughlin phase 
d/ls G [0.9, 1.4] as shown in Fig. [5th), we find a ~ 2.9 ±0.1. Therefore, the universal tunneling 
exponent can be expected if the Laughlin phase without reconstructed edges is confirmed in 
graphene. 

4. Conclusions 

We compare the FQH states at v = 1/3 in models for GaAs-based and graphene-based devices. 
Although the CLL theory predicts a tunneling current I oc V a with a universal a = 3, the 
exponent can be affected by edge reconstruction, which exists in modulation-doped GaAs-based 
samples. However, in a graphene-based device, edge reconstruction can be avoided by using a 
very thin dielectric between the graphene layer and the metallic gate. We demonstrated the 
universal exponent in the edge Green's function calculation. The exponent is always about 3 
once the ground state is in Laughlin region and smaller than 3 outside of the Laughlin region 
which consists with the CLL theory and the experiments in GaAs. 
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